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Abstract 

It is demonstrated that a structured material formed by nonconnected crossed metallic wires may 
enable negative refraction over a wide frequency range. This phenomenon is a consequence of the 
anomalous dispersion characteristics of the material, particularly of the fact that the isofrequency 

O 

contours are hyperbolic. These properties rely on the nonlocal response of the crossed wire mesh, 
and establish a different paradigm for obtaining negative refraction without left-handed materials. 

PACS numbers: 42.70. Qs, 78.20.Ci, 41.20.Jb, 78.66.Sq 
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I. INTRODUCTION 



Negative refraction has undoubtedly been in the spotlight in recent years. Much of the 
interest in this fascinating phenomenon is caused by the fact it contradicts our experience 
that light is bent by common materials (e.g. glass) in a such a way that the projections 
of the incident and transmitted rays along the interface are oriented in the same direction. 
It was first predicted by Veselago [l( forty years ago that materials with simultaneously 
negative permittivity and permeability ("left-handed" materials) may enable negative re- 
fraction. However, such exotic materials were unknown at the time, and only some years 
ago they were made available in the form of metamaterials [2|. Shortly after this finding 
negative refraction was demonstrated at microwaves [3j, but, since structuring materials 
in the nanoscale is still a true challenge, only very recently negative refraction was finally 
revealed at optical frequencies using a truly three-dimensional metamaterial 4|. 

It is well-known that negative refraction may be obtained without left-handed materials. 



For example, one possibility is to use indefinite anisotropic materials, for which one com- 

Another possibility is to engineer the 



ponent of the permittivity is negative jl], 
dispersion characteristic of photonic crystals [10|, lUj] - Here, we demonstrate a distinctively 
different route to obtain all-angle broadband negative refraction using a spatially dispersive 
material formed by a crossed wire mesh. It is shown that due to the strongly nonlocal re- 
sponse of the structured material the dispersion contours of the propagating mode consist 
of two hyperbolas, and that this property makes the group velocity (energy flow) to be 
refracted in an anomalous manner at an interface with air. 

The material considered here is a crossed wire mesh of nonconnected metallic wires with 
radius r w . The orientation of the two arrays of wires is determined by the perpendicular 
unit vectors ui and U2 [Fig. QJi]. It was demonstrated in our previous work 12| that such 
material may be characterized by an anomalously high positive index of refraction in the long 
wavelength limit, and that such property may enable the propagation of very subwavelength 
guided modes, as demonstrated experimentally in [l3]. Recently, we have shown that the 
resonant excitation of such guided modes may enable subwavelength imaging [l^ . 

Evidently, the crossed wire mesh has an anisotropic electromagnetic response. The results 
of our previous studies [3, H, [3] assumed electromagnetic propagation in a plane normal 
(yoz) to the planes of wires. Here, we demonstrate that for propagation along a direction 
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FIG. 1: (Color online) Panel (a): Geometry of the metamaterial: two nonconnected wire meshes 
are oriented at right-angles. The distance (along y) between perpendicular adjacent wires is a/2. 
Panel (b): Isofrequency contours of the fundamental plane wave mode for propagation in the xoz 
plane with electric field in the same plane. The text insets indicate the value of the normalized 
frequency coa/c. The arrows show the orientation of the electric field associated with a given 
isofrequency contour. 

parallel to the planes of wires (xoz) the propagation properties may be dramatically different, 
and may enable negative refraction. To demonstrate these potentials first we will characterize 
the isofrequency contours of the bulk metamaterial. As reported in 
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171 ]. in the 



long wavelength limit the crossed wire mesh may be characterized by the dielectric function 
e = ehU y Uy + £nUiUi + £ 2 2U2U2 (e/i is the relative permittivity of the host medium), with 

En (u, ki) = e h (l- - - J p -V i = l,2lepsWM (1) 



(u/efeu-kfs 

1 /2 

where (3 p = [2n/ (In (a/2nr w ) + 0.5275)] 1 /a, c is the speed of light in vacuum, k = 
(k x , k y , k z ) is the wave vector, ki = k.u,, and a is the distance between adjacent parallel 
wires. For simplicity, it is assumed that the wires are perfect conductors. The effect of 



metallic loss can be easily accounted for using the more general formulae of Ref. 



12j, and 



is relatively small provided the radius of the wires r w is larger than the skin depth of the 



metal 12], which may be verified through the infrared domain. It is convenient to choose 
the system of Cartesian axes so that ui = (1,0,1)^/2 and u 2 = (-1,0, l)/y/2 [see Fig. QJi]. 

Consider the case in which the electromagnetic wave propagates in the xoz plane, k y = 0, 
with magnetic field polarized along y (the time variation e~ tuJt is suppressed): 

H = iZoe^-'Uy , k = (k x , 0, k s ) IH field (2) 

The corresponding electric field is given by: 



Using the Maxwell's Equations it is straightforward to verify that the wave vector must 
follow the dispersion characteristic: 

k 2 k 2 

x — 2 + -2—2 = lldispEq (4) 

k 2 -(u/c) e n k 2 -(u/c) e 2 2 

which, using Eq. epsWM, may be reduced to a polynomial equation of third degree in (3 2 
with (3 = ^/ShUi/c: 

-A(3 6 + 8/3 4 ((3 2 + k 2 x + k 2 z ) - 

1 {Ap 4 p + 8(3 2 k 2 x + 8p 2 k 2 z + 6k 2 x k 2 + 5k x + bk£) + 

(k 2 x - k 2 ) 2 (2(3 2 + k 2 x + k 2 z ) = OldispEQ (5) 

It may be verified that when the operating wavenumber is much smaller than the plasma 
frequency, /3 <C /3 P , there is a unique positive solution for f3 2 . This means that for low 
frequencies and for a given k = (k x , 0, k z ) there is always a propagating mode with electric 
field in the xoz plane. The isofrequency contours of this propagating mode are represented 
in Fig. [lb for a material with r w = 0.05a formed by metallic wires standing in air (e^ = 1). 



As seen in Fig. [Tb and also reported in Ref . 15j, [16] , for a fixed frequency the isofrequency 



contours consist of two hyperbolas with asymptotes running along the directions Ui and 
U2. These hyperbolic contours resemble in part the isofrequency contours of an indefinite 
material, however there is an important difference: the isofrequency contour of an indefinite 



material consists of a single hyperbola 



0. 



In Fig. [Tb we have also represented the electric 



field vector lines, which are qualitatively similar to those in an indefinite material. 
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FIG. 2: (Color online) Angle of transmission of the energy flow (Poynting vector) as a function of 
the angle of incidence, for different frequencies of operation. The radius is r w = 0.05a and the host 
permittivity is £h = 1. The inset represents the geometry of the problem showing the incident, 
reflected and refracted waves. Notice that the wave vector and the Poynting vector are not parallel 
in the metamaterial. 



It is we 



1 known that hyperbolic contours may enable negative refraction at an interface 



with air [5|, [s, Q]. In order to investigate this possibility, we consider the geometry shown 
in the inset of Fig. [2} which shows an incoming plane wave illuminating the metamaterial 
along the direction 9%. The angle of refraction Q t for the energy flow (determined by the 
Poynting vector) can be calculated using the fact that the projection of the wave vector 
onto the interface is preserved [lOj]. Thus, for a given frequency u and angle of incidence 9t, 
the wave vector associated with the transmitted wave is of the form k* = (uj/c sin#j, 0, &*), 
where k\ is calculated using the dispersion equation dispEQ [19|- The Poynting vector is 

n □ 

normal to the associated isofrequency contour [18|, [20(. As in indefinite media, the angle 
between the wave vector and the Poynting vector in the metamaterial is acute. 

In Fig. [2] we show the angle 9t as a function of the angle of incidence 9i, for different 
frequencies of operation. Consistent with hyperbolic shape of the isofrequency contours, the 
angle of transmission 9 t is negative, i.e. the wave group velocity suffers, indeed, negative re- 
fraction at the interface of the crossed wire mesh with air. The results of Fig. [2] indicate that 
this phenomenon (unlike in photonic crystals and left-handed materials) is very broadband, 
being observed for a wide range of frequencies. Notice that at the considered frequencies 
the electrical size of the unity cell of the crossed wire mesh is small, a < A, as required 
in order that homogenization theory can be applied. The emergence of negative refraction 
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FIG. 3: (Color online) Transmission coefficient as a function of normalized k x for a fixed frequency 
and different L. Solid lines: magnitude; Dashed lines: phase. The lattice constant is such that 
uia/c = 0.6, and the radius of the wires is v w — 0.05tt. The curves associated with L — 1.25Aq were 
calculated for a regular dielectric slab with permittivity e = 2. 

in the crossed wire mesh has a simple physical justification. Indeed, consider the scenario 
depicted in FigJIJi, which shows the incoming wave illuminating the metamaterial along a 
direction such that the incoming electric field is nearly parallel to the direction Ui. In these 
circumstances, it is clear from the topology of the metamaterial that the incoming wave will 
interact mainly with the set of wires directed along uj. But since these wires are tilted by 
—45° with respect to the interface, it is apparent that as the wave propagates in the crossed 
wire mesh it suffers a negative spatial shift, or in other words, since the dominant path of 
propagation is expected to be along the wires parallel to ui, the group velocity suffers neg- 
ative refraction. Thus, in a certain sense, each set of wires behaves as a waveguide, and the 
polarization of the incoming wave controls which set of wires is "activated" and which set of 
wires is "dormant" . This heuristic interpretation is of course only a very rough description 
of the complex wave interaction between the crossed wires, but enables one to visualize and 
relate the negative refraction to the microstructure of the material. 

In order to confirm the homogenization results, we used the full wave commercial simu- 
lator CST Microwave Studio™ to calculate the transmission coefficient T of a crossed wire 
mesh slab under plane wave incidence, with magnetic field along the ^-direction [Fig. [3J. 
Specifically, we have calculated T as a function of k x (x- component of the wave vector of the 
incoming wave) for a normalized frequency of operation ua/c = 0.6. Notice that k x is deter- 
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mined by the angle of incidence: k x = u/c sin^. It can be seen in Fig. [3] that the crossed 
wire mesh is relatively transparent to radiation for all incident angles, and that \T\ may be 
close to unity. This property is observed over a very wide frequency range (not reported 
here for brevity); in general, the transmissivity tends to improve with increasing frequency. 
These results indicate that the crossed wire mesh may be well-matched to free-space. 

It is interesting to analyze the variation of the phase of T with k x . For convenience, let 
us write T = |T| e - *^, where <ft — — argT. It can be seen in Fig. [3] that is a decreasing 
function of k x . Interestingly, such behavior is completely different from that in a conventional 
dielectric slab with positive permittivity (green dashed line in Fig. [3]), for which increases 
with k x . Indeed, as proven next, the slope of is equal to the spatial shift A suffered 
by an incoming beam when it crosses the slab [inset of Fig. [3]. In fact, suppose that an 
incoming beam (e.g. with Gaussian profile) impinges on the planar slab. Suppose that the 
magnetic field at the input interface is H y = H l y (x). Using Fourier theory the incoming 
beam can be written as a superposition of plane waves, as H l y (x) = J H y {k x ) e lhxX dk x . 
Evidently, T = T(k x ) may be regarded as the transfer function of the slab. This means 
that the magnetic field at the output is such that: H° (x) = f H y [k x ) T {k x ) e lkxX dk x . But, 
if the spatial spectrum of the incoming beam is highly concentrated at the wave number 
k° x = lo jc sin 6i, i.e. if the beam is a quasi-plane wave propagating along the direction 8i, a 
straightforward analysis (similar to the one used to define group velocity) shows that the 
magnetic field at the output plane is such that: H° (x) ~ T (k®) e lAk °H y (x — A), where 
A = d(f)/dk x . Therefore, apart from a transmission coefficient, the field at the output plane 
differs from the field at the input plane from a spatial shift A, which is completely determined 
by the slope of the phase. This analysis is completely general and is valid for an arbitrary 
material slab. It provides a simple criterion to test the emergence of negative refraction in 
metamaterials, by testing if A is positive or negative. Applying the proposed theory to the 
crossed wire metamaterial, we see from Fig. [3] that the slope of <fi is negative, and thus it 
follows that the spatial shift A is also negative. This is a direct proof of the emergence of 
negative refraction in the metamaterial, and fully supports the homogenization results. In 
particular, for the curve associated with L = 15a (L = 10a) we have numerically calculated 
(for incidence along = 45°), A = d<p/dk x = — 0.5L (— 0.39L), which indicates that the 
transmission angle is 9 t = tan _1 A/L = —27° (9 t = —21°). These values concur well the 
value predicted by homogenization, 6 t = —25° [Fig. [2]. 
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FIG. 4: (Color online) Normalized |E| 2 in the vicinity of the crossed wire mesh calculated by solving 
the Maxwell's equations using the Method of Moments. The incoming wave has a Gaussian profile. 
Both the metamaterial slab and the incoming beam are periodic along y. The wires lie in planes 
parallel to the xoz plane and have radius R = 0.05a. Each plane of wires contains 90 wires. The 
frequency of operation is ioajc = 0.6. Panel (a): L = 10.0a; Panel (b): L = 21.0a. Panel (c): 
L = 31.0a. 

To demonstrate in a conclusive manner the appearance of negative refraction, we used 
the Method of Moments (MoM) to simulate numerically the response of a finite width 
metamaterial slab illuminated by an incoming beam with Gaussian profile. It is assumed 
in the simulation that both the structured material and the incoming beam are periodic 
along the y-direction, with period equal to a. The MoM numerically solves the Maxwell's 
Equations and takes into account all the details of the microstructure of the crossed wire 
mesh, yielding the exact solution of the problem, apart from "numerical noise" . It is assumed 
in the simulations that each plane of inclusions (parallel to the xoz plane) is formed by 90 
wires with radius R = 0.05a. As depicted in Fig. [Tk, the wires in alternate planes are 
perpendicular. The Gaussian beam illuminates the slab along the direction 8i = 45°, and 
has a beam waist equal to 2wq = 4. OA at the normalized frequency of operation ua/c = 0.6 
(i.e. a = 0.1A). The width of the slab along x is approximately W = 90a\^2 ta 12A. In 
Fig. H]we represent the calculated squared amplitude of the electric field (which is roughly 
proportional to the beam intensity) in the vicinity of the crossed wire mesh for different 
thicknesses of the slab. The dashed line represents the propagation path (maximum of the 
electric field amplitude), and clearly shows that the incoming wave is bent with a negative 
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transmission angle at the interfaces with air. The angle 6 t calculated directly from the 
spatial shift A is B t = -22°, B t = -25°, and 6 t = -30° for panels (a), (b), and (c) of Fig. H 
respectively. These values match relatively well the theoretical value 6 t = —25° [Fig. [2]. 

It is interesting to note that the electric field amplitude is significantly lower in the 
metamaterial slab, as compared to the air regions, particularly in panels (b) and (c) of Fig. 
HI One of the reasons is that the field distribution of panels (b) and (c) was calculated at 
a plane equidistant from two adjacent wire planes, whereas the field distribution of panel 
(a) was evaluated on a plane with wires. The second reason is that the wave impedance 

n 

r) in the crossed wire mesh is lower than in free-space j!2lj . and thus, since beam intensity 
is roughly |E| 2 /2?7, the conservation of energy requires that the squared amplitude of the 
electric field is lower in the metamaterial. Despite the difference between the impedance in 
the structured material and free-space, it is evident from the simulations that the level of 
reflections is relatively weak. 

In conclusion, we have demonstrated that a crossed wire mesh may enable negative refrac- 
tion over a wide frequency band. The described phenomenon does not rely on a resonance 
of the inclusions, and due to this reason the effect of loss is expected to be relatively small, 
particularly when the radius of the wires is smaller than the skin depth of metal [12]. A 
block of the considered material enables negative refraction either if the interface is normal 
to x- or to the z- direction. In this regard, the response of the metamaterial is fundamen- 
tally different from that of a conventional indefinite anisotropic material, which only yields 
negative refraction when the interface is normal to the principal axis along which the per- 
mittivity is negative. The described results illustrate the richness of the physics of nonlocal 
materials. This work is supported in part by Fundagao para a Ciencia e a Tecnologia under 
project PDTC/EEA-TEL/71819/2006. 
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